Spin dynamics of the spin-Peierls compound CuGeOs under magnetic field 
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The magnetic field-driven transition in the spin-Peierls system CuGeOa associated with the closing 
of the spin gap is investigated numerically. The field dependence of the spin dynamical structure 
factor (seen by inelastic neutron scattering) and of the momentum dependent static susceptibility 
are calculated. In the dimerized phase (H < H c ), we suggest that the strong field dependence 
of the transverse susceptibility could be experimentally seen from the low temperature spin-echo 
relaxation rate I/T20 or the second moment of the NMR spectrum. Above H c low energy spin 
excitations appear at incommensurate wave vectors where the longitudinal susceptibility Xzz(°) 
peaks. 

PACS numbers: 74.72.-h, 71.27.+a, 71.55.-i 
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Low dimensional quantum magnets have drawn great 
attention in recent years. It is believed that the inorganic 
compound CuGeC>3 is a good realisation of quasi-one di- 
mensional (ID) weakly coupled spin— 1/2 antiferromag- 
netic chains. At low temperature, a Spin-Peierls (SP) 
transition was experimentally established from a rapid 
drop of the spin susceptibilityLH. Simultaneously with 
the opening of a spin gap, the lattice undergoes a struc- 
tural change characterised by a small distortion along 
the chains direction (X-axis) as seen by x-rays diffrac- 
tion studiesB. 

The Peierls instability in S=l/2 Heisenberg chains has, 
on the theoretical side, been a long standing problemtl 
first connected to experiments on organic quasi-lD sys- 
tems such as (TTF)CuS 4 C 4 (CF 3 )4 (Ref. |) where TTF 
stands for tetrathiafulvalene. The new CuGeC>3 com- 
pound is the fist inorganic quasi-lD SP system with a 
critical temperature T$p ~ UK. 

To describe a single Cu 2+ spin- 1/2 chain of CuGeC>3 
we start with the following effective hamiltonianB, 



h = y, J (*) s * • s *+i + J ' E s * • S 



i+2 



(1) 



where i is a site index (i=l,...,L) and Si is an elec- 
tron spin operator. A small dimerisation of the near- 
est neighbor (NN) exchange coupling is due to the lat- 
tice distortion and mimics the effect of the phonons, 

j(i) = j(i+s(-iy). 

Simple quantum chemistry arguments^ suggest that, 
due to side groups and in contrast to the usual case, the 
(almost) 90° Cu-O-Cu superexchange path is antiferro- 
magnetic (i.e. J > 0) with J ~ UOK (11.6 meV). In 
addition, the second nearest neighbor exchange coupling 



J' is also antiferromagnetic and hence leads to frustra- 
tion. A sizeable ratio a — J' / J is suggested by the mag- 
netic susceptibility measurementsEl above T$p (S = 0) 
which can be well fitted by parameters like a = 0.24 and 
J = 150K (Ref. |ora = 0.36 and J = 160K (Ref. |). 

In the dimer (D) or Spin-Peierls phase, for tempera- 
tures below Tsp ~ 14K, the chains become dimerized 
with an alternation of somewhat shorter and longer Cu- 
Cu bonds along the X-axis. Physically, this corresponds 
to the formation of singlet dimers on the stronger bondsfl. 
In the simple (6 — 0) J-J 1 Heisenberg chain the tran- 
sition towards a gapped disordered ground state (GS) 
occurs at a c & 0.24ltL3. The parameters proposed by 
Castilla et al.U hence correspond almost exactly to the 
critical point and, at low temperature, an extra cou- 
pling to the lattice must be advocated to give the correct 
zero tempp;ature finite spin gap seen in neutron inelastic 
scatteringO. This is taken care of by an ad hoc dimeri- 
sation 5 = 0.030. In addition to the value of the spin 
gap such parameters seem also to reproduce satisfacto- 
rilly the whole dynamical spin structure factorEj. How- 
ever, it should be noted that a similar agreement on the 
gap value can also be achieved using the set of parame- 
ters of Ref. H provided a smaller S = 0.014 is usedt 3 ! In 
any case, it is clear that both magnetic frustration and 
spin-phorpji interactions conspire to the formation of the 
SP phasetil. In the following, we shall use both sets of 
parameters. 

The small value of the spin gap (As ~ 2.1meV) 
makes the magnetic field properties of CuGeOa particu- 
larly interesting. The Zeeman energy —gfisSz leads, for 
T < Tsp and for increasing magnetic field (H c ~ 12.5 T 
at low temperature), to a phase transition to an incom- 
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mensurate or solitonic This transition is di- 

rectly connected to closing of the spin gap by the mag- 
netic fieldlia. An accurate theoretical description of the 
full spin dynamics under magnetic field is still missing 
and is necessary for the interpretation of neutron inelas- 
tic scattering experiments under magnetic field or the 
NMR experiments. In the following, we investigate the 
magnetic field behavior of the dynamical spin structure 
factor and of the static q-dependent magnetic susceptibil- 
ity which can be indirectly studied in NMR experiments. 
Dynamical correlations are calculated on finite clusters 
by standard exact diagonalisation methodsta. 




S zz (q,uj,H) = S z 



(2) 



where "to is the S z — ground state and the sum over 
n is, in fact, restricted to the S z = 1 components of the 
excited triplet states. Hereafter, the units of the mag- 
netic field H are such that g\iB = 1< The longitudinal 
part is field independent while the field dependence of 
the tranverse part is simply due to the Zeeman splitting 
of the S=l multiplet producing a shift of the poles of the 
(5-function. 



3 



J 5 





Sz=1 

o=0.36 8=0.014 (a) 


















q= 7i/10 





M/J 




1/2 i 

q 

FIG. 1. (a) Dynamical spectral function S zz (q,uj) calcu- 
lated at momenta q = (L = 28) in the D-phase. From 
bottom to top, n goes from 1 to L/2 (= 14) and the parame- 
ters are indicated on the plot, (b) Momentum-dependence of 
the various structures appearing in S zz (q,u>) for the same pa- 
rameters as in (a). Results for both 20- and 28-sites clusters 
are shown. The open triangles and diamonds are experimen- 
tal data taken from Ref. O. 

In the D phase the transverse and longitudinal dynam- 
ical spin structure factor are given by 

S+_(<Z, u,H) = \(*n\S+\* )\ 2 6(oj — E n + H — E ) , 

n 
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FIG. 2. (a) Dynamical spectral functions S zz (q,ui, H), cal- 
culated at momenta q = 2nJ (L = 20) in the S z — 1 ground 
state (for H > H c ) (b) Momentum-dependence of the two 
main structures with lower energy in S zz (q,u) and upper 
bound of the continuum for the same parameters as in (a). 

Before discussing in more details the role of the mag- 
netic field, we briefly summarize our main results on the 
field independent longitudinal structure factor. The re- 
sults for S zz (q, lj) calculated on a 28-sites chain are shown 
in Fig. 0(a). A sharp low energy feature is seen which has 
been interpreted as a magnon branch or, equivalently as 
a spinon-spinon bound state. A lot of weight is also seen 
above this band. It is interesting to notice that a second 
sharp peak exists above the magnon branch and might a 
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priori be interpreted as a second bound state below the 
continuumEI The excitation energy($f the three lowest 
triplet states as a function of momentum q is shown, for 
two different chain lengths, in Fig. |l|(b) together with the 
upper limit of the continuum of excitations. It is clear 
that, in contrast to the two first triplet states, the exci- 
tation energy of the third one (open symbols) is strongly 
size dependent and, very likely, converges to the value of 
the second excitation energy. This means that the second 
and third excited states belong to a continuum of states. 
We conclude that, for these parameters, there is only a 
single magnon branch separated from the continuum by a 
gap as seen experimentallyEi For completeness, we also 
show a comparison of the magnon dispersion path the 
available experimental data of Regnault et al.tLil. The 
agreement is excellent ,_even slightly better than the fit 
obtained by Haas et al.E3 for the parameters proposed in 
Ref. 0. Note that, although finite size effects are neglige- 
able around q — n/2, a careful finite size scaling of the 
gap at q = ir (see Ref. |l3|) is necessary and, in fact, gives 
a gap very close to the experimental value. 
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FIG. 3. Momentum dependence of static susceptibilities, 
(a) x+- f° r several values of the reduced fields h = H/H c . (b) 
\zz for several values of the magnetisation m z . Parameters 
are indicated on the graph. 



We now turn to the discussion of the effect of the mag- 
netic field H. For H < H c (H c ~ A#), the transverse 

part S-i is simply obtained by a global shift of H of 

the spectral weight of Fig. |I|(a) towards low energies. 
At H c the GS acquires a finite magnetisation S z and 
the lattice modulation becomes incommenturate. In the 
following, we shall use a sinusoidal modulation J(i) = 
J(l + dcos(Ksi)) of momentum K$ = - S z ) 

which minimizes the total energy. However, we stress 
that the effects directly connected, to the change of the 
lattice modulation are very weakE3. In the GS carrying 
a finite magnetisation the longitudinal dynamical spin 
structure factor becomes 

S zz (q,^,H) ^Y,\( xS ^\ q, o(Sz))\ 2 S(io - E n - E ) , (3) 



where the field dependence enters in the GS |*o(5'z))- 
Results are shown in Fig. ||(a-b) for S z = 1. The low en- 
ergy spectral weight appears now away from momentum 
7r at a momentum q m i n (H). Its position (see Fig. §(b)) 
in q-space shifts continuously with the field H according 
to q m in(H) = 7r(l — 2m z ) where m z — S z /L is the mag- 
netization. This is consistent with the mapping of the 
XY model on the non-interacting tight-binding fermion 
model, the change of the magnetic field being related, in 
the fermion picture, to a change of the chemical potential 
and hence of the Fermi momentum kp (Ref. 

Besides inelastic neutron scattering, Nuclear Magnetic 
Resonance (NMRVis a very efficient technique to probe 
the spin dynamicsEj, and in particular, the static (zero 
frequency) susceptibility tensor Xa(q)- The dominant 
component of the interaction between nuclear spins of 
copper in CuGeOa is namely the one mediated by the 
fluctuation of electronic spins, i.e., by the non-local elec- 
tron spin susceptibility. Neglecting the direct nuclear 
dipole - dipole interaction, the second moment (M2) of 
the NMR line-shape obtained-in the standard spin - echo 
experimentEj will in generalE3 depend on both \zz and 

Xh susceptibilities. However, taking advantage of the 

strong anisotropy of the copper nuclear spin hyperfine 
coupling (Azz /Aj_) 2 ~ 100, (where Z is the principal 
axis perpendicular to the d x 2_ y 2 orbital of copper spin),Ej 
we can select only one dominant contribution to M2 by 
applying the magnetic field in the appropriate direction. 
For H || Z or H \\ to the chain axis, we thus measure 
iW^zz or M2,h , respectively, where 



M 2 . a = 



C a {jY, A zzxl(q)~[jY, A2 zzX a (q)} 2 } , 



(4) 



with a — zz or H — , C zz = 0.69(ri7„/f/^ B ) 4 /8?i 2 , C + _ = 
pC zz for the (-1/2, 1/2) transition and C+_ = 9pC Z z/16 
for (±3/2, ±1/2) transitions, where p < 1 is the propor- 
tion of " like^, nuclear spins participating in spin-exchange 
interactionc3. The g-dependence, Azz = Azz{q), of the 
Azz coupling was found to be essentially negligible (< 
10%) with A zz {q = tt) ~ -440 kOeEa 
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The static susceptibilities can be easily obtained from 
the knowledge of the dynamical structure factors, 



X a (q,H)= f 
Jo 



S a (q,uj,H) 



duj . 



(5) 



Results obtained on a 20-sites cluster are shown in Fig. 
||(a-b). In the D-phase (h < 1, m z = 0) Xzz(q) and 

{q) are peaked at momentum q = it. For increasing 

field, Xzz{q) remains constant (x Z z(q,H) = \x+-{q,H = 
0)) while a strong singularity developps at q = n in 

Xh (q, H) when h — * 1. Above the critical point (h > 1), 

the longitudinal part Xzz{q,H) becomes field dependent 
with a maximum at the momentum g m ; n corresponding 
to the appearance of spectral weight at low energies in 
the dynamical spin structure factor. 
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FIG. 4. Second moment Ma,+- vs reduced field H/H c . To 
avoid spurious divergences when h —* 1, the momentum q = n 
is not included in the corresponding sum of Eq. ^j. 

The field dependence of the susceptibilities are ex- 
pected to be directly reflected in the copper NMR data 
in CuGeC>3. The M 2jZZ value is thus predicted to be field 
independent, while M2.^ should strongly grow as we in- 
crease h — ► 1. (The increase of M^h (h) is plotted in 

Fig. U in reduced units, i.e., with C a = 1 = Azz(q), as 
the value of p is not known). Indeed, NMR spectra in the 
D-phase, i.e., the width of NMR lines (at 4.2 K) reported 
in Ref. ^6|, seem to confirm this behavior. These results 
have to be corroborated by further systematic investiga- 
tions. 

In conclusion, we have shown that the strong mag- 
netic field dependence of the spin fluctuations could lead 
to new experimentally observable effects. In particular, 
it is suggested that the behavior of the transverse spin 
susceptibility can be experimentally extracted from the 
second moment of the NMR spectrum. 
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